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Abstract 



> 

Tjlj- ' We compute the single-particle inverse lifetime, along with the conductivity- 
(N ■ 

' derived scattering rate, for a metallic system in an s-wave superconducting 



state. When both electron-phonon and electron-impurity scattering are in- 



' eluded, we find that while these scattering rates are in qualitative agreement, 

a" 

in general quantitative agreement is lacking. We also derive results for the 



quasiparticle lifetime within the BCS framework with impurity scattering, 
which makes it clear that impurity scattering is suppressed for electrons near 



^ ' the Fermi surface in the superconducting state. 

I. INTRODUCTION 

The quasiparticle lifetime is a concept which is useful for clarifying the nature of the 
system of interest, i.e. is it a Fermi vs. marginal Fermi vs. a Luttinger liquid, etc (Bedell 
et al, 1990). Various techniques are available for measuring lifetimes (Kaplan et al. 1976); 
perhaps the most direct is through tunnel junction detection (Narayanamuti et al. 1978). In 
this paper we analyze electron scattering rates, as measured by microwave and far-infrared 
conductivity measurements, and examine their relationship to quasiparticle lifetimes. We 
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will reserve the name "quasiparticle lifetime" to refer to the single-particle property to be 
defined technically in the next section. In contrast we will use the term "scattering rate" or 
"scattering lifetime" to refer to a property derived from a response function, such as the op- 
tical conductivity in this case. The experimental results for YBaCu307_^(Bonn et al. 1992, 
1993) and their theoretical implications (Bonn et al. 1993; Berlinsky et al. 1993; Klein 1994; 
Bonn et al. 1994) have been discussed extensively in the literature. In these works the 
microwave conductivity was used to extract information about quasiparticle scattering; it is 
commonly and often tacitly presumed that the "lifetime" (inverse scattering rate) derived 
from these sorts of measurements is closely related to the quasiparticle lifetime, as derived 
from the single-particle Green function. In this work we explore this relationship and com- 
pute both the single-particle lifetime and the two-particle scattering rate, including both 
electron-impurity scattering and electron-phonon scattering. The latter involves inelastic 
scattering processes and introduces significant complications. We use an Eliashberg formal- 
ism to analyze this problem, with some further approximations which will be made clear 
in the relevant sections. A similar approach has already been used to analyze the normal 
state problem (Shulga et al. 1991; Marsigho and Carbotte 1995). Here, certain aspects 
of quasiparticle lifetimes are clarified and the analysis is extended to the superconducting 
state, following our recent work on the optical conductivity (Marsiglio and Carbotte 1995; 
Marsigho et al. 1996). 

We will proceed by carefully defining the single-particle lifetime. We review two ways in 
which conductivity data can be used to extract a scattering rate, one based on the two-fluid 
model, and the other based on a straightforward Drude fit to the low frequency conductivity. 
The instances in which these procedures yield a qualitative or quantitative facsimile of the 
quasiparticle inverse lifetime will be clarified. In this way, one can evaluate the usefulness of 
the two-fluid hypothesis, for example, in systems with both elastic and inelastic scattering 
channels. 

This paper is divided into two sections. The first reviews quasiparticle lifetimes, which are 
calculated from the single-particle Green function. In the normal state the single-particle 
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equations for electron-impurity scattering (in the Born approximation) and electron-phonon 
scattering are well known (Mahan 1981; Allen and Mitrovic 1982; Grimvall 1981). In the 
superconducting state, we first examine the BCS case with electron-impurity scattering. 
By "BCS" we mean the limit where the pairing interaction is instantaneous, and in this 
case, local, i.e. the gap function is independent of frequency and momentum. The single- 
particle spectral function is given analytically, showing explicitly the role of the so-called 
BCS coherence factors which are somewhat disguised in the clean limit. When impurity 
scattering is included, such an expression should be used in lieu of the more commonly 
adopted phenomenological form found in standard texts (Schrieffer 1983). When inelastic 
electron-phonon scattering is included in the problem, perturbative methods within the 
Eliashberg framework are used. This follows closely the work in Kaplan et al. (1976); 
however, here we discuss the quasiparticle lifetimes in the presence of both impurity and 
phonon scattering. 

The second section elucidates two methods recently used to extract scattering rates from 
optical and microwave conductivity measurements, and how these relate to the quasiparticle 
inverse lifetime. We should make it clear at the start that the calculation we use for the 
optical conductivity omits vertex corrections. The main result of this is that the "1 — cos ^" 
factor that occurs in a Boltzmann formulation of transport properties is absent (Mahan 
1981), so in fact, when theoretical comparisons between the quasiparticle lifetime and the 
conductivity-derived scattering rate are made, the agreement will in general be better than 
it would be, given an exact calculation. Conversely, poor agreement between these two 
quantities in the theory would almost certainly result in poor agreement between the two 
measured quantities. 

II. QUASIPARTICLE LIFETIMES 

Quasiparticle lifetimes in clean electron-phonon superconductors were discussed long ago 
by Kaplan et al. (1976). Before proceeding to their relationship with the optical conductivity 
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we give a brief review of the formalism (Kaplan et al. 1976; Nicol 1991). In the normal 
state the quasiparticle has energy and lifetime defined by the pole of the single particle 
retarded Green function, i.e. the zero of 

G'^{k,uj + i5) = uj - ek-T.{uj + i5), (1) 

where S(u; + i5) is the electron sclf-cncrgy. More specifically the solution ujk, where uj^ = 
Ek + ^{oJk + defines the quasiparticle energy E^, and inverse lifetime by 

uJk = Ek- iTk/2 (2) 

Similarly, in the superconducting state, the inverse quasiparticle lifetime is defined by (twice) 
the imaginary part of the pole in the single-particle Green function. The diagonal component 
of the single-particle Green function is (Scalapino 1969) 

u^Z^{uj) -ei- 02 (w) 

where Z{ijj) and (t){uj) are the renormalization and pairing functions given by solutions to 
the Eliashberg equations (Eliashberg 1960; Scalapino 1969; Marsigho et al. 1988) which are 
repeated here for convenience: 

(t>m 



(f){(jj) = ttT ^ [A(a; - iUm) - ij,*{uJc)9{uJc - \uJm\)]- 



^Juj'^{ljJ — u) — (j)'^{uj — u) 
+[7V(.) + /(. + ^)]^=^^^^—=\ (4) 



and 



UJ[UJ) — U + lirl \[UJ — lUJm 



+m dv a^F{v) J [N{v) + f{v - u)] 



Ujl^Z^{iUm) + (p m 



-y/cl;2(u; — I/) — (p'^^uj — v) 



+ ^<- + -' . (5) 
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where oj{uj) = ujZiuS). Here, iV(z/) and f{y) are the Bose and Fermi distribution functions, 
respectively. The electron-phonon spectral function is given by o?F{v\ its Hilbert transform 
is A(z). The Coulomb repulsion parameter is fi*{uJc) with cutoff Uc- A negative value for 
this parameter can be used to model some BCS attraction of unspecified origin. The renor- 
malization and pairing functions are first obtained on the imaginary axis at the Matsubara 
frequencies, i.e. u = iun = mT{2n — 1), with 0^ = (piiujm) by setting the complex variable 
u in these equations to the Matsubara frequencies (Owen and Scalapino 1971; Rainer and 
Bergmann 1974). Then the equations are iterated as written, with uo set to a frequency on 
the real axis. Note that the square roots with complex arguments are defined to have a 
positive imaginary part. 

Actually, these functions are obtained at frequencies just above the real axis, and then, in 
principle, the pole is given by the zero of the denominator continued to the lower half-plane. 
However, since the solutions are readily known only along certain lines in the complex plane 
(e.g. the imaginary axis or just above the real axis) we follow previous authors and look for 
the pole perturbatively. That is, we write uo = E — iV, and linearize the imaginary part so 
that (Scalapino 1969; Kaplan et al. 1976) 

^^^^ = ^^(Er ~ Ezm 

Note that E is determined by equating the real parts (again linearizing in imaginary com- 
ponents), which yields 



E 



\ zKE) ■ ^'^ 

In the normal state, 0(c<j) = 0, and Eq. (^ simplifies to the known normal state result 
(Allen and Mitrovic 1982; Grimvall 1981). The quasiparticle energy and lifetime are given 
by Eqs.(^ and (^, which yield in the low energy limit (Allen and Mitrovic 1982; Grimvall 
1981): 



T{E = 0) 



i + 47r dua"^ 



F{u)[N{u) + f{u)] 



(9) 



1 + A*(T) 



where /(z^) is the Fermi function and 



oo 



:2F(z/)/mV^'(- + ^^) 



(10) 



Typical resuhs for various electron-phonon spectral functions are shown in Fig. 1 for the 
clean limit over a wide range of temperatures. 

Returning to the superconducting state, at the Fermi surface efc = so Eq. (|^) gives 
E = (f)i{E)/Zi{E) = Ai{E), which becomes the definition for the lowest energy excitation, 
i.e. the gap in the excitation spectrum. It has become common practice to consider E as an 
independent variable, and then to study T{E) as a function of E. In Fig. 2 we show T{E) vs. 
E for various temperatures, using a Debye model spectrum. It is clear that the scattering 
rate near the gap edge (shown by the arrow) decreases very quickly as the temperature 
decreases. Note that r(_E') actually becomes negative at intermediate temperatures near 
E ^ 10 meV. This is real (i.e. not a numerical artifact) and is simply a property of the 
function T{E) given by Eq. (|^) through the linearization procedure. The true pole in Eq. 
(0) will always have a negative imaginary part (i.e. F is positive). 

What happens when both the BCS limit and the clean limit are taken ? Then Z2{E) — > 
and 02 (-E') 0, i.e. the functions involved in the solution are pure real, with Zi{E) 1 and 
(f)i{E) — s> A(T), where A(T) is obtained self-consistently from the BCS equation. Clearly 



then the quasiparticle energy is = yef, + and the scattering rate, F = 0. Thus, within 

the BCS approximation in the clean limit the quasiparticle states are infinitely long-lived, 
as there is no means by which a quasiparticle can decay. With either impurity or phonon 
scattering, quasiparticle decay becomes possible. In this way the BCS approximation in the 
clean limit is pathological. 

Let us first examine the situation where only electron-impurity scattering is present, i.e. 
no phonon scattering occurs. Then the pairing and renormalization function become 
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Z{u) = Z^iico) + ^^^^=— (11) 

where the subscript 'cl' refers to the clean hmit and 1/r is the (normal) impurity scattering 
rate. The gap function is defined 

A(cu) = (l){u)/Z{u) (12) 

and is independent of impurity scattering. In the BCS limit, 

i Asgncj 



(t){L0) = A + 



2r ^u'^ - A2 



Z{u) = l + ^-^^. (13) 
As before, using the perturbative approach, the quasiparticle energy is given by = 



e| + A2, but the scattering rate is now (Pethick and Pines 1986) 



1 _ A2 1 



r = -1^—-^ = - ^ . (14) 

r \E\ r^^. 

Eq. ( P^ ) shows that at the Fermi surface impurities are ineffectual for electron scattering, 
i.e. this is another manifestation of Anderson's theorem (Anderson 1959). This happens 
abruptly at Tc; as soon as a gap develops the scattering rate becomes zero. Away from the 
Fermi surface the scattering rate is reduced from its value in the normal state. In particular, 
well away from the Fermi surface the scattering rate in the superconducting state approaches 
the normal state rate. This behaviour is illustrated in Fig. 3. 

The calculations in Fig. 3 rely on a perturbative search for the pole in the lower half- 
plane, as given by Eq. (^). A more rigorous calculation reveals poles "within the gap", but 
these do not appear in the spectral function because of coherence factors. This can be seen 
by rewriting Eq. (|^) as 

^ ,1 \ 1 A UJ \ 1 

Gu{k,uj)= - 1 + 



2 V Vc^2 _ ^2 j ^^2 _ ^2 _ + ^sgna; 
i 1- r ] , — (15) 



The coherence factors, [1 ± ^^al^2 j ; such that the spectral function, A(k, cu) = 
— ^ImG{k,uj), is zero for \uj\ < A, independent of the impurity scattering rate, as can 
be readily verified explicitly from Eq. (|1^). On the other hand the precise pole of the Green 
function is given by solving for uj = E ~ iT /2 in the two coupled equations 

- (r/2)2 = + 4 - (16) 
r = ^; (17) 



For 1/r = the solution is as before, E = = ^ e\ + , and F = 0. For finite impurity 
scattering, however, the solution depends on at the Fermi surface (e^ = 0): 



E = Ja^-{^)\ F = for ^<1 (18) 



E = 0, F = ^(i)2-(2A)2 for ^>1 (19) 

In either case the solution has a real part that lies within the gap, that is, the quasiparticle 
residue at the pole is zero. The "relevant" energies (i.e. for which the residue is non-zero) are 
\E\ > A, and then the scattering rate given by Eq. (|17[) agrees with the linearized solution, 
Eq. (^. However, the quantity T{E) is only physically meaningful when it is evaluated at 
the pole energy, which occurs below the gap value, so in the "relevant" region T{E) is merely 
a well-defined function. Away from the Fermi surface the quasiparticle energy requires the 
solution of a quadratic equation and E may lie below or above the gap, depending on and 
1/r. Once more solutions with \E\ < A have zero spectral weight. In Fig. 4a (b) we show 
the real (imaginary) parts of the pole for some typical parameters. As can be seen, some 
atypical behaviour can occur as a function of temperature, for example when = 0. As 
the temperature is lowered the pole moves from a point on the negative imaginary axis to 



the origin (near T/Tc = 0.8) and then moves along the real axis towards E = \/A^ — 



' J_\2 
^2r/ ■ 



Nonetheless, comparison of Fig. 4(b) with Fig. 3 shows that the perturbative is in qualitative 
agreement with the non-perturbative solution. 
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Coherence factors are known to play an important role for two-particle response functions 
(such as the NMR relaxation rate or the microwave conductivity). In particular they lead to 
the Hebel-Slichter (Hebel and Slichter 1959) singularity in the NMR relaxation rate. However 
they have largely been ignored in single-particle functions. Eq. ( p^Sf ) shows, however, that 
the coherence factors play a very important role in the single-particle spectral function, in 
that they maintain a gap equal to A, even when impurities cause the single-particle pole to 
have a real part whose value lies in the gap. 

The frequency dependence of the spectral function is shown in Fig. 5. For energies close 
to the Fermi surface the spectral function is dominated by the square root singularity at the 
gap edge, which comes from the coherence factors rather than the single-particle pole. For 
larger energies the peak present is due to the pole and differs from the normal state spectral 
function only at low frequencies. 

We return now to the strong coupling case. In the superconducting state, we do not 
have access to the analytic continuation of the Green function to the lower half-plane. (It is 
true that we can attempt a representation through Fade approximations, as has often been 
done in the past, to analytically continue functions from the imaginary axis to the real axis. 
However, these are notorious for producing spurious poles (Marsiglio et al. 1988).) We thus 
follow Karakozov et al. (1975), and use an expansion near = 0. We also follow Kaplan et 
al. (1976) and henceforth use Eq. (|^) for the scattering rate, noting that while the energy 
E at which T{E) is evaluated ought, in principle, to be computed self-consistently as was 
done in BCS, i.e. eqs. ( p!6| , p!7D , in practice we will choose the relevant energy, and so treat 
energy as an independent variable. 

Karakozov et al. (1975) pointed out that at any finite temperature the solutions to the 
Eliashberg equations have the following low frequency behaviour: 

Z{uj) ^ + ^ (20) 
A(cj) ^ diuj"^ - i52UJ (21) 

where Zi,72,5i and 82 are real (positive) constants. This implies that a quasiparticle pole 
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exists with 



r ^ ^(4vr r dva'F{y)g{v)[N{v) + /(z/)] + ^) 
Zi ^ JO r ^ 



(22) 



(23) 



where we have used Eq. (^) and assumed — > (i.e. we are near the Fermi surface). In 
Eq. ( p3D (^(z/) is the single electron density of states in the superconducting state, 



which is non-zero for z/ = at any finite temperature (Karakozov et al. 1975; Marsiglio and 
Carbotte 1991; Allen and Rainer 1991). In Fig. 6 we show the quasiparticle scattering rate, 
r{T,E = 0) vs. T/Tc for various impurity scattering rates, 1/r, in the superconducting 
state. The normal state result is also shown for reference. Note that in the clean limit there 
is an enhancement just below in the superconducting state, but at low temperatures 
there is an exponential suppression, compared to the power law behaviour observed in the 
normal state. When impurity scattering is present there is an immediate suppression below 
Tc (a knee is still present, however, in the superconducting state). It is clear that at low 
temperatures the superconducting state is impervious to impurity scattering, as one would 
expect. 

Is r(T, E = 0) the relevant inverse lifetime by which properties of the superconducting 
state can be understood ? Strictly speaking the answer is no, particularly at low tempera- 
tures, where the spectral function essentially develops a gap for low energies. In Fig. 7 we 
plot the spectral function at the Fermi surface, A^kp, uj) for various temperatures (Marsiglio 
and Carbotte 1991). While no true gap exists at any finite temperature, it is clear that for 
low temperatures the spectral weight at low frequency is exponentially suppressed. Thus 
we have a situation similar to that in BCS theory, where the quasiparticle pole occurs at 
an energy where the spectral weight is essentially zero. It is more relevant to inquire about 
the scattering rate (as defined by Eq. (|])) evaluated above the "gap-edge" where a large 




(24) 
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spectral weight is present, and quasiparticles are more likely to be populated (as in detector 
applications). When impurities are added the spectral peak is broadened, even at low tem- 
peratures, as shown in Fig. 8 for an intermediate temperature, T/Tc = 0.5. However the 
lineshape is very asymmetric as a gap remains at low frequencies (particularly prominent at 
low temperatures). 

Fig. 7 demonstrates that the relevant energy is a function of temperature. In Fig. 9 we 
plot r{T,E = A(T)) vs. T/Tc, where A(T) is determined from the relation (at the Fermi 
surface) 

E = ReA{E,T). (25) 

By Eq. (pSf ) we understand that the E = solution is excluded (similarly the very low 
energy solution is also excluded). We are interested in the conventional solution which gives 
rise to the peak in the spectral function illustrated in Fig. 7 or 8. If no non-zero solution 
is present (as is the case near T^) then we utilize the E = solution. Note that in this 
case care must be taken when obtaining the E ^ limit of Eq. (^. Also there is present 
in this definition a discontinuity at some temperature near T^, which is where a nonzero 
solution first appears. In this way we hope to show the scattering rate at an energy where 
the spectral weight is large, and therefore of most relevance to observables. At any rate, it is 
clear by comparing Fig. 9 to Fig. 6 (see also Fig. 2) that there is very little difference in the 
scattering rate in the gap region of energy. However, as the energy increases beyond the gap 
the scattering rate increases, resulting in short "lifetimes" , even at zero temperature (Kaplan 
et al. 1976). Thus, while we have argued that the quasiparticle lifetime corresponding to 
the pole of the Green function is of limited use because the quasiparticle residue there is 
zero, in practice the close correspondence between Fig. 9 and Fig. 6 illustrates that it 
remains a useful indicator of the scattering rate for energies up to about the gap energy. 



III. THE OPTICAL CONDUCTIVITY 
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A. Extraction of a Scattering Rate from the Conductivity 



Using the Kubo formalism (Mahan 1981), the optical conductivity can be related to a 
current-current correlation function. The final result for the frequency dependence of the 
conductivity in the long wavelength limit is (Nam 1967; Lee et al. 1989; Bickers et al. 1990; 
Marsiglio et al. 1992) 

a(z.) = l(n(z/ + ^5) + — ), (26) 

where the paramagnetic response function, n(i/ + i5), is given in the previous article (Mar- 
siglio and Carbotte 1997), and will not be repeated here. 

To summarize the theoretical framework within which we are working: for a given model 
for the spectral density, a^Fi^u) and a choice of impurity scattering rate 1/r, we can compute 
the conductivity cr(z/) at any frequency and temperature, including effects due to both elastic 
and inelastic scattering mechanisms, the latter being determined by the choice of electron- 
phonon spectral density. As explained in the introduction, however, vertex corrections are 
omitted. 

In the analysis of experimental data it is possible to use several methods to extract 
a single temperature dependent scattering rate. One such method utilized the microwave 
conductivity in YBaCu307_x(Bonn et al. 1992, 1993) and in Nb (Klein 1994), and adopted a 
two-fluid model description of the superconducting state. It was assumed that the absorptive 
component of the conductivity (real part of a at finite frequency, denoted by cri(z/)) was due 
only to the normal component of the fiuid. Thus an expression of the form 

r(r) 



ai(i/, T) 



m m*{T) 



1 



1 + {ut{T)Y 



(27) 



A2(T)_ 

is assumed to hold approximately. In Eq. (^) r(T) has units of a scattering time and A(T) is 
the penetration depth at temperature T. Here we will calculate o"i(z/, T) for a model a^F{v) 
and 1/r using the full expression (pU]). At the same time we can calculate the penetration 
depth, either from a zero frequency limit of Eq. (^), or directly from the imaginary axis 
(Nam 1967; Marsiglio et al. 1990): 
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The idea is to examine the zero frequency hmit of Eq. (^7|), and thus define a scattering 
rate relative to the rate at Tc (Klein 1994): 

r{T) \ \\T)) a,{T) ' ^ ^ 



where it has been assumed that the mass enhancement factor in Eq. (|27|) does not change 
with temperature. In this way an effective scattering rate was extracted from microwave 
data for YBaCu307_^(Bonn et al. 1992, 1993) and Nb (Klein 1994). 

To see how closely this scattering rate follows the quasiparticle scattering rate, we show 
results in Fig. 10 for t(Tc)/t(T) defined by Eq. ( p9D in the superconducting state (solid 
curve) with which we compare the inverse quasiparticle lifetime for both the normal (short- 
dashed curve, Eq. ), and the superconducting (long-dashed curve, Eq. (||) ) states. The 
results are based on a Debye model spectrum for q;^-F(z/) used in the previous section, with 
mass enhancement parameter A = 1 and Tc = 100 K. (A negative fi* is required.) Note that 
the normal state scattering rate (short-dashed curve) is only sensitive to the low frequency 
part of a^F(z/) at low temperatures: a u'^ dependence in a'^F{v) implies a dependence 
in l/r(T) (Again, vertex corrections would alter this to the familiar law (Ashcroft and 
Mermin 1976).) The agreement between the inverse lifetime (long-dashed curve) and the 
scattering rate defined by the two-fluid model (solid curve) in the superconducting state is 
remarkable, although later we shall see that this quantitative agreement occurs only with 
this particular Debye spectrum. This indicates that the two-fluid description makes sense 
(Berlinsky et al. 1993; Bonn et al. 1994) at least qualitatively, and the quantities shown in 
Fig. 10 apply to the normal component of the superfiuid. Fig. 10 illustrates the comparison 
in the clean limit, where the two-fiuid description is expected to be most accurate (Bonn et 
al. 1994). Before investigating impurity dependence, we turn to a second possible procedure 
for extracting a scattering rate from conductivity data (Romero et al. 1992; Tanner and 
Timusk 1992), which is simply a generalization of that used by Shulga et al. (1991) to the 

13 



superconducting state. One simply fits the low frequency absorptive part of the conductivity 
to a Drude form: 

X ne^ 1 'r*(T) , , 

m m* /m 1 + [yT*y 

As described in Shulga et al. (1991) and Marsiglio and Carbotte (1995), it is possible to fit a 
Drude form to the low frequency part of the optical conductivity in the normal state. Such a 
fit is also possible in the superconducting state (Romero et al. 1992). Theoretically, the fit is 
problematic in a BCS approach because there is a Hebel-Slichter logarithmic singularity at 
low frequency at all temperatures in the superconducting state. However, with the Eliashberg 
approach, the Hebel-Slichter singularity is smeared, and one can fit a Drude form over a 
limited range of frequency. Such a fit for a Debye spectrum is also included in Fig. 10 (dotted 
curve). While the fit is in qualitative agreement with the inverse lifetime, this method of 
characterizing the inverse lifetime in the superconducting state is clearly not as accurate as 
the two-fluid model. The fits themselves are shown in Fig. 11. It is clear that the fits fail 
at sufficiently high frequency, as one would expect, but that they characterize well the low 
frequency response in the superconducting state. 

Fig. 10 clearly shows that there is considerable freedom and hence ambiguity in extract- 
ing a temperature dependent scattering rate from conductivity data. Nevertheless, in the 
clean limit it is evident that the two-fluid model is a useful device for extracting the low 
energy inverse quasiparticle lifetime in the superconducting state. 

With the addition of impurities the situation changes considerably. This is illustrated 
in Fig. 12, where the same calculations as in Fig. 10 are shown, but with an additional 
impurity scattering, 1/r = 25 meV, included. The use of formula (0), inspired by the 
two-fluid model, gives a scattering rate (solid curve) that falls much less rapidly around 
T = Tc then does the inverse quasiparticle lifetime in the superconducting state (long- 
dashed curve). The latter curve drops almost vertically as the temperature drops below Tc, 
as has already been discussed. The result based on the Drude fit (dotted curve) shows a 
peak which is reduced in size from that in the clean limit (Fig. 10). Indeed, for increased 
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impurity scattering, the peak just below Tc disappears. In both cases a rapid suppression is 
expected just below T^. in the case of the inverse lifetime, this suppression is a consequence 
of Anderson's theorem (Anderson 1959), as already discussed. In the case of the Drude fit, 
it is easy to see that this is the case in the dirty limit. In the dirty limit the conductivity 
is almost flat as a function of frequency on the scale of l/r, at T^. Just below however, 
a gap in the spectrum begins to develop, so that weight is shifted from roughly the gap 
region to the delta function at the origin. Any low energy flt will then use a Lorentzian 
width which senses this depression in the conductivity, which is on an energy scale of the 
gap. This represents a significant suppression from the normal state scattering rate (infinite 
in the dirty limit). Here we are in an intermediate regime, with l/r = 25 meV (note: 
A(r = 0) = 20.2 meV). The corresponding fits for Fig. 12 are shown in Fig. 13. 



B. Possible Origin of the Low Temperature Conductivity Peak 

As an aside we wish to further examine the low frequency conductivity vs. reduced 
temperature, a quantity measured in the high Tc oxides by microwave techniques (Nuss et 
al. 1991; Bonn et al. 1992, 1994). In these experiments a peak was observed in the real 
part of the low frequency conductivity as a function of temperature, the origin of which has 
been discussed in many contexts. Here, we add yet another possibility. In Fig. 14 we show 
the real part of the conductivity, cri(t/) vs. reduced temperature, for several low frequencies. 
We continue to use a model Debye spectrum to simply estabhsh qualitative effects. Note 
the relative insensitivity to frequency, a feature of strong coupling pointed out in Marsiglio 
(1991). Also note the lack of a coherence peak just below Tc- Nonetheless, a broad peak 
exists at lower temperatures, somewhat reminiscent of that observed in YBaCu307_j;(Nuss 
et al. 1991; Bonn et al. 1992, 1994). This peak exists because of a competition between an 
increasing scattering time (making a increase) and a decreasing normal component (making 
(7 decrease, particularly as T — > 0). We should note that this peak is most prominent in 
the clean hmit. As Fig. 13 indicates (see values at the intercept), the peak is absent for a 
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sufficiently large impurity scattering rate. 

It is of interest to examine what dependence these results have on the electron-phonon 
spectral function. As an extreme we utilize a spectrum which is sharply peaked at some high 
frequency, and, in contrast to the Debye model employed above, coupling to low frequency 
modes is absent; such a spectrum models a strong coupling to an optic mode. We choose a 
triangular shape for convenience, starting at luq = 34.8 meV with a cutoff at lje = 35.5 meV. 
The coupling constant is chosen so that the mass enhancement value is A = 1, in agreement 
with that chosen for the Debye spectrum. As was the case there, a negative /i* is used to 
give Tc = 100 K, and the zero temperature gap was found to be close to the Debye value. 

In Fig. 15 we plot the real part of the low frequency conductivity, cri(i/) vs. reduced 
temperature, now using the triangular spectrum for a^F{v). In contrast to the result for 
the Debye model, the low frequency conductivity is strongly frequency dependent at low 
temperatures. In fact, for i/ = 0, the conductivity appears to diverge. We believe that 
at sufficiently low temperature, this curve will actually achieve a maximum and approach 
zero at zero temperature, but we have been unable to obtain this result numerically. Once 
again there is a competition between an increasing scattering time and a decreasing normal 
density component as the temperature is lowered from T^. Here, however, the increasing 
scattering time appears to be overwhelming the decrease in normal fluid density. The key 
difference with the Debye spectrum is that here the spectrum has a big gap, so that the 
lifetime is increasing exponentially with decreasing temperature, already in the normal state. 
Recall that in the Debye case the increase followed a power law behaviour with decreasing 
temperature. Since the decrease in normal fluid density is always exponential, this term 
dominates in the Debye case, whereas in the case of the gapped spectrum the competition is 
subtle, and will depend strongly on the details of the spectrum (an electron-phonon spectrum 
with a much smaller gap will yield a zero frequency conductivity which approaches zero at 
zero temperature, for example). 

The physics of this conductivity peak is different from what has been already proposed 
for YBaCuaOy-^. It has been suggested that the excitation spectrum becomes gapped due to 
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the superconductivity, i.e. a feedback mechanism exists which creates a low frequency gap in 
the excitation spectrum as the superconducting order parameter opens up below T^. Such a 
scenario has been explored within a marginal Fermi liquid scheme (Nuss et al. 1991; Nicol 
and Carbotte 1991a, 1991b), and is also consistent with thermal conductivity experiments 
(Yu et al. 1992). Here the conductivity peak arises because the spectrum is already gapped, 
and the scattering rate is sufficiently high at Tc because Tc itself is fairly high. We should 
warn the reader that this mechanism requires a "fine tuning" of the spectrum, i.e. a large 
gap is required in the phonon spectrum. An alternate boson spectrum may be operative, 
and then a less extreme boson spectrum may work, if the normal state electron density falls 
at a correspondingly slower rate as the temperature is lowered. This would be achieved 
if the superconducting order parameter has nodes, for example (Schachinger and Carbotte 
1996). 

Note that for any non-zero frequency the conductivity has a visible maximum, and 
quickly approaches zero at sufficiently low temperature. Nonetheless this turn around occurs 
for yet another reason: as one lowers the temperature the Drude-like peak at low frequencies 
gets narrower while at the same time the magnitude of the zero temperature intercept 
increases. For any given finite frequency, then, a temperature is eventually reached below 
which this frequency is now on the tail of the Drude-like peak. This means that while the 
zero frequency conductivity increases, that at any finite frequency will eventually decrease 
as the width becomes smaller than the frequency. So in this case the increase in scattering 
lifetime still dominates the decrease in normal fluid density, but, because we are fixed at a 
finite frequency, the conductivity decreases. Note that frequencies of order 0.01 meV 2.4 
GHz) are within the range of microwave frequencies that are used in experiments. 

To show this more explicitly we illustrate in Fig. 16 the various scattering rates obtained 
with the gapped spectrum. These are to be compared with those shown for the Debye model 
in Fig. 10. Clearly the quantitative agreement between the scattering rate inspired by the 
two-fluid model and the inverse lifetime that was seen in Fig. 10 with the Debyc spectrum 
was fortuitous. While a qualitative correspondence between these two entities continues to 
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exist with the gapped spectrum, they are no longer in quantitative agreement. We have 
verified that this is generically true, by investigating other spectra, not shown here. 

IV. SUMMARY 

We have investigated the quasiparticle lifetime in an Eliashberg s-wave superconductor, 
generahzing earlier work (Kaplan et al. 1976) to include impurity scattering as well. We 
find that the quasiparticle lifetime becomes infinite at low temperatures, independent of 
the impurity scattering rate, which we understand as simply a manifestation of Anderson's 
theorem (Anderson 1959). Thus, on general grounds, within a BCS framework, the inverse 
quasiparticle lifetime should collapse to zero in the superconducting state. We have also 
emphasized the importance of the so-called coherence factors for the single-particle spectral 
function, when impurity scattering is included within the BCS case. 

Two methods have been investigated for extracting the scattering rate from the low fre- 
quency conductivity. One relies on a two-fluid model picture, and the other simply utilizes a 
low frequency Drude fit. Neither should necessarily correspond very closely to the quasipar- 
ticle inverse lifetime, and we find that in general they do not, quantitatively. Qualitatively, 
however, the scattering rate defined by either procedure gives the correct temperature de- 
pendence for the inverse lifetime. In the presence of impurities, the two-fluid prescription 
appears to be less accurate, presumably because such a prescription takes into account only 
the lower normal fluid density as the temperature is lowered, and not the fact that impurity 
scattering is less effective in the superconducting state. Thus we caution that the inter- 
pretation of a conductivity- derived scattering rate as a quasiparticle inverse lifetime, while 
qualitatively correct, is quantitatively inaccurate. 
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FIGURES 
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FIG. 1. (a) Mass enhancement parameter, as defined by Eq. (|lO| ) in the text, vs. reduced 
temperature T/Tc, for various electron-phonon spectral functions. In all cases Tc = 100 K. The 4 
spectra used are a Debye spectrum (solid line), linear spectrum (dotted line), a spectrum propor- 
tional to ^/u (dashed line), and a triangular spectrum (dot-dashed line). In all cases the strength 
is such that A = 1. In the first three cases a cutoff frequency equal to 30 meV was used. In the 
last case, the spectrum starts at 34.8 meV and is cut off at 35.5 meV. (b) Inverse lifetime, T(T) (in 
meV) vs. reduced temperature for the same spectra as in (a). Note the different low temperature 
behaviour, depending on the low frequency characteristics of the electron-phonon spectrum used. 
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FIG. 2. Scattering rate, T{T,E) (in meV) vs. E (in meV) for various temperatures in the 
superconducting state. The zero temperature gap at the Fermi surface is indicated by the arrow. 
Note that for T/Tc = 0.5 the function plotted actually becomes negative (near 10 meV). The 
physical pole occurs at an energy, E, where r(r, E) is always positive, however. 



26 



1.0 - 



0.5 



0.0 



1 1 1 

1/(tA) = 1 


1 


e/A = 5 

/ 

t 




/ 

_ 0.5 




0.1 


0.0 


1 1 1 



0.0 



0.5 

T/T 



1.0 



FIG. 3. Scattering rate normalized to the zero temperature gap, r(T, £')/A vs. reduced tem- 
perature T/Tc, for various quasiparticle energies, e. These are computed using the perturbative 
expansion, Eq. (jlj). The impurity scattering rate is l/(rA) = 1. Note that on the Fermi surface 
(soHd curve), the scattering rate is zero immediately below Tc- 
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FIG. 4. (a) Real and (b) imaginary parts of the quasiparticle pole vs. reduced temperature, 
T/Tc, with 1/(tA) = 1 and various quasiparticle energies. These are computed non-perturbatively 
from Eqs. ([T6| , |T7|) . Note that at the Fermi surface the quasiparticle energy has an abrupt onset at 
a temperature somewhat below Tc. The results in (b) are similar to those in Fig. 3, except that at 
the Fermi surface, the decrease in scattering rate below Tc is not as abrupt. 



28 




29 




30 



FIG. 5. The single-particle spectral function, A{k,uj), vs. normalized frequency, to/ A, with 
1/(tA) = 1, for various quasiparticle energies, e/A. The Fermi surface result (solid curve) has 
particle-hole symmetry. The result for e/A = 0.5 would have a peak within the gap (between -1 
and 1) except that the coherence factors in Eq. ( p!5| ) give zero residue for the gap region. As the 
quasiparticle energy increases, the spectral function begins to resemble the normal state spectral 
function. Small square-root singularities still exist, nonetheless at the particle and hole gap edges. 
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FIG. 6. The quasiparticle scattering rate, r(r, = 0) vs. T/Tc for various impurity scattering 
rates, 1/r, in both the superconducting and normal (long-dashed curves) states. In the clean limit 
(solid curve) there is an enhancement immediately below T^. When impurity scattering is present, 
this scattering is immediately suppressed in the superconducting state, as shown by the dotted 
and dashed curves. Below a temperature of about O.STc the scattering rate is independent of the 
amount of impurity scattering present in the normal state. A Debye electron-phonon spectrum 
was used with A = 1 and cutoff frequency, a;£) = 30 meV. 
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FIG. 7. The single-particle spectral function, A{kF, Oij) at the Fermi surface, vs. frequency, in 
the clean limit, for various reduced temperatures. The spectral function is considerably broadened 
near Tc, due to temperature alone. At the two lowest temperatures shown, while there is no true 
gap in the excitation spectrum, this plot makes it clear that, practically speaking, an effective gap 
in the excitation spectrum is present. A Debye electron-phonon spectrum was used with A = 1 
and cutoff frequency, = 30 meV. 
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FIG. 8. The single-particle spectral function, ^(^^,0;) at the Fermi surface, vs. frequency, for 
various impurity scattering rates, as indicated. Note the broadening which occurs with increasing 
impurity scattering. However, spectral weight remains absent in the "gap region" . Results are for 
a temperature T/Tc = 0.5, and with a Debye electron-phonon spectrum with A = 1 and cutoff 
frequency, a;£) = 30 meV. 
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FIG. 9. The quasiparticle scattering rate, T{T,E = A{E)), evaluated at the quasiparticle 
energy, given on the Fermi surface by E = A{E), vs. T/Tc for various impurity scattering rates, 
1/r, in both the superconducting and normal (long-dashed curves) states. These results are in 
quantitative agreement with those in Fig. 6, since the energy scale, A, is still small compared 
to other (phonon) energy scales in the problem. In particular, below about O.STc, the scattering 
rate is independent of the amount of impurity scattering present in the normal state. A Debye 
electron-phonon spectrum was used with A = 1 and cutoff frequency, a;^) = 30 meV. 
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FIG. 10. Various normalized scattering rates vs. reduced temperature. The normal and su- 
perconducting scattering rates come from the quasiparticle inverse lifetime, given by Eq. @, at 
zero energy. The two-fluid result comes from Eq. ( p9|) while the Drude fit is obtained by fitting 
Eq. (^) to the low frequency conductivity in the superconducting state. Note the agreement of 
the scattering rate as extracted from the two-fluid analysis with the inverse lifetime in the super- 
conducting state. A Debye electron-phonon spectrum was used with A = 1 and cutoff frequency, 
a;£) = 30 meV. 
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FIG. 11. The low frequency conductivity in the superconducting state (sohd curves) along with 
their fits based on Eq. (|30|) (dashed curves). These fits were used in Fig. 10 (dotted curves). 
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FIG. 12. Same as Fig. 10, except now with an impurity scattering rate of 25 meV. Note that 
the two-fluid analysis agrees poorly with the quasiparticle inverse lifetime. 
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FIG. 13. The low frequency conductivity fits used in Fig. 12. 
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FIG. 14. The very low frequency conductivity as a function of reduced temperature, in the clean 
limit. A Debye electron-phonon spectrum was used with A = 1 and cutoff frequency, w^i = 30 meV. 
Note that the results are relatively insensitive to frequency (the v = 0.01 meV result, given by the 
dotted curve, is essentially hidden by the zero frequency result). A microwave experiment yields 
essentially zero frequency results. 
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FIG. 15. Same as for Fig. 14, but now with the triangular spectrum as described in the text. 
Note that the zero frequency conductivity appears to diverge as T — >■ 0. A microwave experiment 
will yield a large low temperature peak as a function of reduced temperature, whose magnitude will 
depend strongly on the frequency. These results are for the clean limit. The peak is also reduced 
as impurity scattering is added (not shown). 
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FIG. 16. Comparison of scattering rates (as in Fig. 10) calculated for the triangular spectrum, 
in the clean limit. Note that the normal state result approaches T = exponentially due to the 
gap in the a^F{i') spectrum (in Fig. 10 the corresponding curve approached zero with a power 
law behaviour). While the results in the superconducting state are qualitatively similar, they no 
longer agree quantitatively with one another, as in Fig. 10. 
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